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The intermediate and late-time behaviour of a massive Dirac field in the background of static 
spherically symmetric dilaton black hole solutions is investigated. The intermediate asymptotic 
behaviour of a massive Dirac field depends on the mass parameter as well as the wave number of 
the mode, while the late-time behaviour has a power law decay rate independent of both. 

PACS numbers: 04.50.-|-h. 



I. INTRODUCTION 

The late-time behaviour of various fields in the spacetime of a collapsing body is of a great importance in black 
hole physics. Regardless of details of the collapse or the structure and properties of the collapsing body the resultant 
black hole can be described by just a few parameters such as mass, charge and angular momentum, black holes have 
no hair. The manner and rate with which the hair of the black hole decays is thus an important question. In what 
follows we begin by reviewing some of the old an new work on this problem. 

Price in Ref.[l| for the first time studied the neutral external perturbations. He found that the late-time behaviour 
is dominated by the factor ^-(2^+3)^ fgj. each multipole moment I. The decay along null infinity and along the 
future event horizon the fall off was found '^l to be like that and where u and v were the outgoing 

Eddington-Finkelstein and ingoing Eddington-Finkelstein coordinates. In Ref.[3| the scalar perturbations on the 
Reissner-Nordstrom background for the case when | Q \< M was studied. A late time dependence like ^-(2^+2) -^^as 
found, while for | Q \= M the late-time behaviour at fixed r is governed by t~^^^^\ Charged scalar hair decayed 
slower than a neutral one while the late-time tails in gravitational collapse of fields in the background of 

Schwarzschild solution was reported by Burko Q- The intermediate and late-time pattern of hair decay was also 
considered in Reissner-Nordstrom background in Ref.jl]. The very late-time tails of massive scalar fields in the 
Schwarzschild and nearly extremal Reissner-Nordstr/"om black holes were the subject of 0, It was found that 
the oscillatory tail of a scalar field decays like t~^^^ at late times. Power-law tails in the evolution of a charged massless 
scalar field around the fixed background of a dilaton black hole were studied in [ll[ , while the case of a massive scalar 
field was treated numerically in Ref. ■ The analytical proof of the decay pattern both for intermediate and late-time 
behaviour was presented in [3] ■ In Ref. [ij] the late-time tails of massive scalar field were studied in the spacetime of 
stationary axisymmetric black hole and it was found that the power law index of —5/6 depended neither on multiple 
mode I nor on the spin rate of the considered black hole. All the above cases involved bosonic fields. 

On the fermionic side, the problem of the late-time behaviour of massive Dirac fields were studied in the spacetime 
of Schwarzschild black hole ^15^, while in the spacetime of an Reisner-Nordstrom black hole was analyzed in [TB|. 
The case of the intermediate and the asymptotic behaviour of charged massive Dirac fields in the background of 
Kerr- Newman black hole was elaborated in Ref. ■ 

Interest in unification schemes such as superstring/M-theory has triggered an interests in the decay of the hair of 
n-dimensional black holes. As far as the n-dimensional static black holes is concerned, the no-hair theorem for them 
is quite well established [l8| . The decay mechanism for massless scalar hair around higher n-dimensional Schwarzshild 
spacetime was worked out in p^ . The late-time tails of massive scalar fields in the spacetime of an n-dimensional 
static charged black hole was treated in [20| and it was revealed that the intermediate asymptotic behaviour of the 
considered field had the form t~('+"/2-i/2) , The above pattern of decay was confirmed numerically for the case of 
n — 5 and n = 6. One should also mention the results of Ref.[2l|, where the authors obtained fermion quasi normal 
modes for massless Dirac fermion in the background of higher dimensional Schwarzschild black hole. Recently, there 
has been also some efforts to study the late-time behaviour of massive scalar fields in the background of black holes on 
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brane [2^ having in mind the idea that our universe is only a submanilfold on which the standard model is confined 
to, inside a higher dimensional spacetime (23| . On the other hand, the massless fermion excitations on a tensional 
3-brane embedded in six-dimensional spacetime were studied in [23 |. 

The main purpose of the present paper is to extend our knowledge of the behaviour of fcrmionic fields and to clarify 
what kind of mass-induced behaviours play the dominant role in the asymptotic late-time tails as a result of decaying 
the massive Dirac hair in the background of a four-dimensional dilaton black hole. That is in a spherically symmetric 
solution of the low-energy string theory with arbitrary coupling constant a. 

The paper is organized as follows. In Sec. II we gave some general remarks concerning behaviour of Dirac spinors in 
a curved background. Sec. Ill will be devoted to the analytical studies of intermediate and late-time pattern of decay 
of the hair in question, while in Sec. IV we conclude our investigations. 



II. THE DIRAC EQUATION IN A CURVED SPACETIME 

Because of the full spectrum of neutrinos, and their masses and mixing properties is not known, it seems worth 
while summarizing the general situation. Especially as there has recently been some controversy [l^ [l^, [H, about 
the so-called Majorana and Dirac masses and their consequences in gravitational field. 

In four spacetime dimensions we can always use a representation in which the gamma matrices are real and we can 
take the components of all classical fermion fields to take values in a Grassmann algebra over the reals. 

The most general Lagrangian for k four-component Majorana fermion i — 1,2, ... k is thus 

^^JW^Ji^' - V^^MyV^" , (1) 

where ^ = and is the covariant derivative = 9^ + ji^'^^'^^a^b, M and a are tangent and spacetime indices. 

There are related by e°, a basis of orthonormal one- forms. The quantity w"'' = w"'' are the associated connection 
one-forms satisfying de° + Wfc" A e** = 0. On the other hand, 7^ are Dirac matrices satisfying {7*^,7^} — 277"''. We 
are using a mainly plus metric signature convention. In four spacetime dimensions, for example, the gamma matrices 
may be taken to be real. The 4k x 4fc matrices T = Tij + ^^Oij, Mij = fiij + ^^Vij with 75 = 70717273, 75 = —1 and 
Tij,i2ij,Vij symmetric and atj anti-symmetric matrices which we take to be independent of time and position. 
There is an action on the spinors by GL{k,C) which satisfies 

^'^S'j^P\ (2) 

where 5* j = exp(a* j + 75/3* j) and we are thinking of 75 as a complex structure on M*'' = C^*^. Elements of C^*^ are 
Weyl (or chiral) spinors for which 75 = z. Under this action the result yields 

T^S^fTS, M^S'MS. (3) 

Using the freedom ^ one may set T — diag ± 1, ±1, . . . ±. So as to have positive energy we demand that all signs 
are positive. Now S G U{k,C) and we may choose it to make M diagonal with real non-negative entries jJMil- Thus 
we arrive at k uncoupled Dirac equations of the form 

j^S/^-m^ij^O. (4) 

If one iterates the Dirac equation and uses the cyclic Bianchi identity in a curved space one gets the following: 

- V^V' + ii^V' + = 0. (5) 

There is no gyro-magnetic coupling between the spin and the Ricci or Riemann tensors [Hj . We see in these calculations 
no sign of the effect claimed in Refs. p6l [27l. [2^ . This is consistent with the equivalence principle, according to which 
all particles should fall in the same way in a gravitational field. Of course, if the matrices T and M were depended 
upon position then, things could be different. In the presence of a dilaton and axion field, this might happen. In 
this paper we shall just consider mass terms. Although the detailed calculations above assume that spacetime is 
four-dimensional, they are readily extended to higher dimensional spacetimes. 



III. THE DECAY OF DIRAC HAIR IN THE BACKGROUND OF A BLACK HOLE SOLUTION 

The treatment of fermions in spherically symmetric backgrounds may be greatly simplified by recalling a few basic 
properties of the Dirac equation. These allow a rapid reduction of the problem to the behaviour of a suitable second 
order radial equation. We shall begin by giving a discussion valid for all spacetime dimensions n (when S0{3) is 
replaced by SO{n — 1)) but our detailed decay results will apply only to the case n = 4. 
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A. Some useful properties of the Dirac operator 

As we saw above, we may assume that the massive Dirac equation in a background metric is given by Eq.(j4])- The 
basic properties of the Dirac operator = 7^V^ on an n-dimensional manifold that we shaU need are 

• for a metric product 

gfj.^dx'^dx'' = gab{x)dx''dx^ + g„iniy)dy"^dy" , (6) 

it decomposes as a direct sum 

P=P.+Py (7) 

• Under a Weyl conformal rcscahng given by 

g^Lv = ^^gfii^, (8) 

it follows directly that we have 

= i: = n'^^"~-^^. (9) 

For a conformo-static metric of the form 

ds^ = -A^dt^ + ^^dx'dx' , (10) 
where A = A{x'^) and $ = $(a;'), i = 1, 2 ... n — 1, we write 

ds^ = A^ (^-dt^ + {^^dx'dx^^ , (11) 

and consequently find the following: 

= A-5("+2)('/at+^)v;, (12) 



where ^ is the Dirac operator of the metric {^^^dx'^dx^ and ijj = A^^" ^^i/j. Now we use the conformal property 
again. One obtains the relation 

/ 4 \ 

i(«-2). 



(13) 



with ^ = Uj v^. 



Since every spherically symmetric metric is conformally flat, a special case of the theory above, is a static metric 
of the form given by 

ds^ ^ -A^dt^ + B^dr"^ + C^dEl_2, (14) 

where A = A(r),B — B(r),C — C{r) are functions only of the radial variable r, and the transverse metric dYj\_2 
depends neither on t nor on r. 

Let us suppose now, that 5* is a spinor eigenfunction on the [n — 2)-dimensional transverse manifold E. Namely, ^' 
satisfies the relation of the form as 

p^-^ = (15) 
By virtue of the properties given above one may assume that the following is satisfied: 

pip — mil), (16) 
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and set what follows: 



^-ir^x«*. (17) 



It can be verified by the direct calculations that Eq. p^ provides the result as 



where we have denoted 



(7"9t + 7'5,)x = A(m- A)^. (18) 
B 

dy = jdr, (19) 

the radial optical distance (i.e., the Regge- Wheeler radial coordinate). Gamma matrices 7^,7^ satisfy the Clifford 
algebra in two spacetime dimensions. 

An identical result may be obtained if a Yang-Mills gauge field is present on the transverse manifold S, but 
now one gets 

I^SA * = (20) 
where p-^ ^ is the Dirac operator twisted by the the connection A^^ . 
Assuming that ip oc e~"^*, one obtains the second order equation for x 



^+uj\ = A\m--)\. (21) 
In what follows, the detailed form of the spinor harmonics and the eigenvalues will not be important. 



B. The Background 

In four spacetime dimensions, the action for the dilaton gravity with arbitrary coupling constant implies 



R - 2V^0Vp 



(22) 



where cj) is the dilaton field, a coupling constant while i^^^ — 2V[^Ai,] is the strength of U{1) gauge field. 
The static spherically symmetric solution of the equations of motion are given by the following line element: 



ds^ = - 1 - ^ 1 - — 



r_ \ 1+ 



dr'^ 



1-^1(1-^ 



R\r)dn^, 



(23) 



where R^{r) = yl — ^ j , while r+ and r_ are related to the mass M and the electric charge Q of the black 
hole 



-2a<l 



r 



2M = r_ 



1 + a- 



1 + q;2 



(24) 



The metric is asymptotically flat in the sense that the spacetime contains an initial data set {Y,end,gij,Kij) with 
gauge fields such that Se„d is diffeomorphic to R'^ minus a ball and the following asymptotic conditions are fulfilled: 



l5y - 5,J\+r\^ag^^ + ■■■^-A^a,...a,g^^+r\K,,\ + ... + r''\da,...a,K,,\ < ^ , 

\F<,p\ + r\daF^p\ + ... + r''\da,...a,F^i3\ < O 



o 



(25) 
(26) 
(27) 



where Kij is the exterior curvature, (f>o is a constant value of the scalar field. 
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C. Spinor No-Hair theorems 

The properties of static spinor fields around Schwarzschild and Kerr black holes and the consequent no-hair prop- 
erties have been investigated by many people including [H, 113, [H, . The basic idea is to study solutions of the 
static Dirac on the background of the black hole. One either considers the case when there are no fermionic sources 
outside the horizon or one constructs a Green function. In the massless static spherical case it is clear from our work 
above that this is equivalent to solving the flat space Dirac equation [s^] 

yS.T/. = 0. (28) 
This may have regular solutions on the horizon. On the other hand, in our case, we have the following: 

On the horizon A = and, unless the solution is extreme, $ ^ 1371 . The extreme case is exceptional because ^ = 
at the horizon in such a way that the spinor ■0 remains finite [37|. Of course, in the non-extreme case, one should 
check that some scalar spinorial invariant blows up, but this can easily be done. 

D. Decay of Fermionic Hair 

We shall now analyze the time evolution of massive Dirac spinor field in the background of dilaton black hole by 
means of the spectral decomposition method. In it was shown that the asymptotic tail is connected with the 

existence of a branch cut situated along the interval —m < lu < m. An oscillatory inverse power-law behaviour of 
massive Dirac field arises from the integral of Green function G{y,y']Lu) around branch cut. The time evolution of 
massive Dirac field may be written in the following form: 



x{y,t)= / dy' 



G{y,y';t)xt{y',0) + Gt{y,y';t)x{y',0) 



(30) 



for t > 0, where the Green's function G{y, y'; t) is given by the relation 

Q2 



, V 



G{y,y';t)^S{t)6(y~y'). (31) 



In what follows, our main task will be to find the dilaton black hole Green function. Using the Fourier transform |25| 
G{y,y';uj) = dt G{y,y']t)e'^^* one can reduce equation (pij) to an ordinary differential equation. The Fourier's 
transform is well defined for Im lo while the corresponding inverse transform yields 

G(2/, y'; t) ^ — dco G{y, y'; iv)e-'^\ (32) 

J -oc+ie 

for some positive number e. By virtue of the above the Fourier's component of the Green's function G{y,y']uj) can 
be written in terms of two linearly independent solutions for homogeneous equation. Namely, one has 



d^ 
dy^ 



+ LJ^-V]x^ = 0, *-l,2, (33) 



where V = A? \ m — ^ 
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The boundary conditions for Xi described by purely ingoing waves crossing the outer horizon of the static 
charged black hole xi — e^*'^^ as y — )■ — oo. On the other hand, X2 should be damped exponentially at z+, namely 
X2 ^ e"^""-'^"?' at y oo. 

Let us assume that the observer and the initial data are situated far away from the considered black hole. In order 
to rewrite Eq. (|33p in a more convenient form we change variables 

X. = ^ (34) 



1-^ 1-^ 



2(l + Q^) 
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where i = 1,2. Then, one can expand Eg. ((33)) as a power series of r±/r neglecting terms of order 0{{uj/r)'^) and 
higher. Under this assumption we reach to the following: 

2cLi^(r_|_ + air^) — m?{rj^ + Q!ir_) + 2\m{l + r_|_) 



2 2 I 



(35) 



— 2Amr_(ai + a2) + m^air_|_r_ 



where we have denoted ai 
It can be verified that Eq, 



and a2 = ^ 



may be solved in terms of Whittaker's functions. Consequently, two basic solutions 
are needed to construct the Green function, with the condition that | a; |> m. Namely, the Whittaker's functions 
Xi — Ms_fi{2Cjr) and X2 — Ws^ji{2Cor) have the following parameters: 



uj'^{r+ + air_) + Am(l + r+) - ^{r+ 



(36) 



2 2 
TO — CiJ . 



On the other hand, the spectral Green function takes the form as 



Gc(r,r';t) = ^ 
zvr 

_ 1 

2^ 



dw 



Xi(r,a;e") X2(?-', '^e'"*) Xi(r,a;) X2(?-','i) 



l^(t;)e'^«) 



l^((i) 



(37) 



where W{Cj) is the Wronskian. 

Further, we focus our attention on the intermediate asymptotic decay of the massive Dirac hair, i.e., in the range of 
parameters M <^ r ^ t <^ M / [mMY . The intermediate asymptotic contribution to the Green function integral gives 
the frequency equal io Co — 0{^ra/t), which in turns implies that (5 <C 1. Having in mind that 5 results from the 1/r 
term in the massive scalar field equation of motion, it depicts the effect of backscattering off the spacetimc curvature 
and in the case under consideration the backscattering is negligible. Taking into account all the above and the fact 
that Cor -^1 and M{a, 6, z) = 1 as z tends to zero, we obtain the resulting expression for spectral Green function 



G (r r'-t) 23^'"^ r(^2A) r(i + A)r(A + 1) 
' ' ^ mV^ Mr(2M)r(i-/i) 

In the limit when t ^ 1/m it implies 

23M-ir(-2A) r(i + A)r(/i + i) 



„(2A+l)7ri 



(rr')^ 



Km 



Ji 



Gcir,r';t) 



Ar(2A) r(i - a) 



1 + e(2'i+i)'^^ ) (rr')5+'^ to'^ r^-^" cos{mt ~ -{jl + 1)). 



(38) 



(39) 



Eg. ((39)) depicts the oscillatory inverse power-law behaviour. We remark that in our case the intermediate times of 
the power-law tail depends only on p, which in turn is a function of the multiple number of the wave modes. 

The different pattern of decay is expected when k 3> 1, for the late-time behaviour, when the backscattering off 
the curvature is important. Consequently, f{oj) when k 3> 1 may be rewritten in the following form: 



r(l + 2/i) r(l-2/l) 
~ 2/i 

(r(i + 2M))2 r(-2/2) r(i + /2-^) 
2^r(2A) n^-p-s) 



(40) 



J2;l(\/8^) j2f,{Vm^' 



where we have used the limit Af5^^(2a)r) « r(l + 2fx){2Cjr)^ 5^i^ J^{\/85ujr). The first part of the above Eg. ((40)) the 
late time tail is proportional to t^^ and it occurs that we shall concentrate on the second term of the right-hand side 
of Eq. (HP)) . It turned out that for the case when k ^ 1 it may be rewritten in the form as 



Gc (2){r,r']t) 



^ / dw e'(2-^— *) e^^, 



(41) 
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where we have used the following definition: 



On the other hand, M yields 

(r(i + 2m))2 r(-2M) 



M 



2fi r(2/i) 



{rr'y 



(42) 



(43) 



At very late time both terms e™* and e^'^^ are rapidly oscillating. From this fact it follows directly that the spinor 
waves are mixed states consisting of the states with multipole phases backscattered by spacetime curvature, which 
most of them cancel with each others which have the inverse phase. Thus, one can find the value of Gc (2) by means 
of the saddle point method. The saddle point integration allows us to evaluate the accurate value of the asymptotic 
behaviour. Namely, it could be found that the value 2^5 — wt is stationary at the value of w equal to the following: 



TT (cj^ (r+ + Q!ir_) + Am(l + air_) — ^-^{r+ + air_) 



(44) 



Evaluating Eq. (PT|) by means of the saddle point integration we achieve finally to the form of the spectral Green 
function for massive Dirac spinor hair. It implies 



Ge(r,r';t) 



^mV3 (,)l 



2m^(r++air_)+2ATO(l+r_|_)— m^(r++air_) 



[mt) 6 sin[mt) x{r,m) xi'i'' T'fn). (45) 



The above equation provides the main result of our calculations. It illustrates the fact that the late-time asymptotic 
decay pattern of massive Dirac hair in the background of spherically symmetric dilaton black hole is proportional to 
-5/6. 



IV. CONCLUSIONS 



In this paper we have treated the problem of the asymptotic tail behaviour of a free Dirac field in the spacetime of a 
spherically symmetric charged black hole solution of dilaton gravity with arbitrary coupling constant a. This theory is 
related to the the low-energy limit of the heterotic string theory and on its own is a generalization of electromagnetism 
by adding scalar field dilaton and coupling constant between t/(l) gauge field and scalar field. 

The resultant intermediate asymptotic behaviour depends on the field parameter mass as well as the wave number 
of the mode. But this is not the final pattern of decay of the massive Dirac hair. Resonance backscattering off the 
spacetime curvature dominates at late times. We have calculated analytically that the pattern of decay in question 
is proportional to t~6. The same result one gets studying the late-time behaviour of free massive scalar fields in 
the same background. One should remark that the above considerations are also applicable to the case of extremal 
dilaton black hole, i.e., to the case when = r_. Thus having in mind Eq. ([l5l) one gets the exact form of the spectral 
Green function for the late-time behaviour of massive Dirac hair for the extremal dilaton black hole in the theory 
with arbitrary coupling constant a. 
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